Can inflationary models of cosmic perturbations evade the secondary oscillation test? 
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We consider the consequences of an observed Cosmic Microwave Background (CMB) tempera- 
ture anisotropy spectrum containing no secondary oscillations. While such a spectrum is generally 
considered to be a robust signature of active structure formation, we show that such a spectrum can 
be produced by (very unusual) inflationary models or other passive evolution models. However, we 
show that for all these passive models the characteristic oscillations would show up in other observ- 
able spectra. Our work shows that when CMB polarization and matter power spectra are taken into 
account secondary oscillations are indeed a signature of even these very exotic passive models. We 
construct a measure of the observability of secondary oscillations in a given experiment, and show 
that even with foregrounds both the MAP and PLANCK satellites should be able to distinguish 
between models with and without oscillations. Thus we conclude that inflationary and other passive 
models can not evade the secondary oscillation test. 
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I. INTRODUCTION 

VO ' 

qq . Secondary oscillations in the angular power spectrum of CMB anisotropies are broadly understood to be a signature 
of inflationary models for the origin of the primordial perturbations. The inflationary models are part of a wider class 
^\ , of 'passive' models, where the perturbations are set up at very early times, and then evolve according to linear 
>> \ ' perturbation theory for many orders of magnitude of expansion. 

In linear perturbation theory, the physics of gravitational collapse (on scales greater than the Jeans length) results 
O ,■ in a 'squeezing' effect as the presence of one growing and one decaying solution draws any initial conditions into a 
prefered region of phase space. On scales below the Jeans length there is oscillatory behavior at early times, but in 
passive models each mode has already experience a long period of squeezing by the time the oscillatory behavior kicks 
in. The region of phase space prefered by the squeezing corresponds to a fixed phase of oscillation in the oscillatory 
epoch. The fact that different wavenumbers have their oscillations fixed at different phases ultimately results in the 
kJ ' characteristic 'secondary oscillations' in the CMB angular power spectrum 
.£h ! The active models are fundamentally different B. The canonical examples are the cosmic defect models. The key 
difference is that these models have a matter component (e.g. cosmic defects) which is evolving in a non-linear manner 
from very early times. One can think of the non-linear component as a source which constantly seeds pertubations 
in the other matter. This kind of non-linear 'force term' disrupts the squeezing effects which are present at the 
linear level, and those realistic defect models for which calculations have been performed have shown no secondary 
oscillations whatsoever in the CMB anisotropy power H-fjJ . 

Thus the detection, or otherwise, of these acoustic oscillations is cited as a strong indication of whether or not the 
evolution of the perturbations was active or passive. In fact, the absence of these oscillations is regarded as one of 
the few ways that aspects of inflationary cosmology can be falsified || . But what if inflation produced a primordial 
spectrum with oscillations, such that the CMB temperature power spectrum has no oscillations? Can a Ci spectrum 
with no oscillations really rule out inflation? In this paper we investigate this possibility. 

Suppose that the temperature spectrum measured by PLANCK turns out not to have oscillations. Note that as we 
see in Section |V], the new Boomerang || and Maxima Jig ] data can be fit equally well by spectra with and without 
oscillations, so we have no indication either way as yetjJ^With future observations we can hope to decide whether 



1 A different issue is whether the single peak that is clearly evident in the current data is too sharp to be produced by defect 
models. One of us answers this question in the affirmative in fill]. 
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or not the temperature spectrum has oscillations. What could we say about an inflationary model that produced a 
temperature spectrum without oscillations, from our other observations, those of the polarization and matter power 
spectra? We look at a model which has oscillations in the primordial power spectrum producing a temperature 
anisotropy spectrum with no oscillations, to sec how this affects the other spectra. We argue that such an inflationary 
model is forced to have polarization and matter spectra which are easily distinguishable from those of active models, 
at least when the best data sets are finally in. Thus, even for the special borderline cases we consider here the presence 
or absence of secondary oscillations remains a robust test of the fundamental nature of the primordial perturbations. 



II. OSCILLATIONS IN PRESENT DAY POWER SPECTRA 



In this section we review the origin of the oscillations in the angular spectrum for passive models, and discuss 
what factors can affect them. We develop our discussion based on the approximate scheme of Hu and Sugiyama [QJ|] 
which is more intuitive and gives a good enough account of the the physics for our immediate purposes. For the full 
calculations which are reported later in the paper we use the full Boltzmann code CMBFAST [ 12| . 

The temperature anisotropy At = AT/T is expanded in plane waves and Legendre polynomials [ |l3| , |l~| . 

A T (x,n,r) = y"d 3 ke lkx A T (k,n,r) 

= /rf 3 ke lk ' x ^(2Z + l)(- l )'A Ti (k,T)P ; (k-n). (1) 
^ i 

The Boltzmann equation must be solved for the photon phase space distribution to get Ar/(k, t). 

The observed sky is expanded in spherical harmonics A^(x = 0, n, r ) = ^2i m a (m^/m( n ); so to compare the theo- 
retical predictions to observations it is convenient to calculate the angular power spectrum, 

Cti = (\aim\ 2 ) 

= (4ir) 2 ( k 2 dkPi(k)\A Tl (k,T )\ 2 (2) 



where < AT/(k, r)A Tl (k', t) >= S(k — k')Pj(fe)|A<n(fe,r)| 2 , P,(fe) is the primordial spectrum, Ap;(fc,T) is the photon 
transfer function and To is the conformal time today. Homogeneity and isotropy mean that these quantities depend 
only on the modulus of the wave vector k. It is important to note that the photon transfer function does not depend 
on the initial power spectrum. 

The angular anisotropy spectra for polarization and the cross-correlation between temperature and polarization are 
defined similarly, 



Ccx = (4tt) 2 W J k 2 dkP t {k)A Tl {k lT() )A El {k,T ). (3) 

The AEi(k, r) are rotationally invariant quantities related to the Api(k,r) that appears in the Boltzmann equation 
by a spin lowering operator As;(fc,r) = — 3 2 Ap/(/c,t) 

The Boltzmann equations in Fourier space for the gauge-invariant scalar temperature and polarization anisotropics 
are|Jl,|l| 

A T + (ikfj, + k)A T = - ik^ + k(A T0 + ifiVb + ^P 2 (n)n), 

A P + {ikfi + k)A P = k Q(l - P 2 (/i))n) . (4) 

where /x = k • n, n is the optical depth, '5 and $ are gravitational potentials, V& is the baryon velocity and n = 
At2 + Apo + Ap2 . Since there is much less power in the polarization anisotropy than in the temperature, II « At2 > 
and can be ignored in the temperature equation as the temperature monopolc and baryon velocity dominate during 
scattering. 
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In order to understand the behaviour of the anisotropics, the visibility function he K can be approximated by a 
delta-function at the time of recombination t* . Solving the Boltzmann equation gives us 0] 

A T/ (fc,r ) a (A T0 + *)(fe,r,)(2I + l)ii(fcfo - r,)) 

+ A T1 (fc,r*)(0'i-i(fc(T - t„)) - (I + l)j;+i(fc(T - t„))) 

+ (21 + 1) J ° dr^ - m(k(T - t)) (5) 

A^(fc,r ) « g(2Z + l)A T2 (fc,n) i^^p (6) 

There are three terms that come from the surface of last scattering. The first two terms contributing to the 
temperature perturbation are the first two multipoles A^o + ^ and Ayr, which come from the intrinsic density 
fluctuations and velocity fluctuations on the last scattering surface respectively. These are created during the tight 
coupling epoch before last scattering, where the photons and baryons are coupled via Compton scattering. The Bessel 
functions act as projectors from fluctuations in k on the last scattering surface to fluctuations in / that we see today. 
The quadrupole At2 which produces photon polarization is created during last scattering, as the electrons and protons 
recombine. 

The third term contributing to the temperature fluctuation is the integrated Sachs- Wolfe effect (ISW), where 
photons travelling towards us are perturbed by fluctuations in the gravitational field. Large angles are dominated 
by the ISW effect, small angles by the contributions made during scattering. Oscillations in the power spectrum are 
created on small scales (inside the horizon at last scattering) as shown below, so the ISW effect will not be relevant 
here. 

With adiabatic initial conditions, on small scales during tight coupling the monopole and dipole terms are [jj] 

A T0 (fc, r) w (1 + R)-V 4 A T0 (k, 0) cos(fcr s (r)), 

A T i(fc,r) a -\/3(l + i?r 3 / 4 A T0 (fc,0)sin(fcr s (T)), (7) 

where R is the ratio of baryons to photons normalised to 3/4 at equality and r s is the sound horizon. We can set 
Aro(fc 7 0) to 1 for all k as the initial conditions are specified in the initial power spectrum. The dipole (velocity 
contribution) oscillates out of phase with the monopole (density contribution), and with reduced amplitude. So the 
temperature anisotropy today, ATi{k,To), is dominated by the density fluctuation at last scattering. 

From the solution to the Boltzmann equation (|5|) we can see that at times shortly after r* the quadrupole behaves 
like At2 ~ An, so the behaviour of the polarization anisotropy today is dominated by the velocity fluctuation during 
last scattering. 

So for adiabatic initial conditions the temperature perturbations on the surface of last scattering are dominated by 
a cosine oscillation, and the polarization perturbations are dominated by a sine oscillation. The bessel functions act 
as projectors from Fourier mode k to multipole I, so for large I we only need to consider large k contributions to the 
angular power spectra to see their oscillatory structure, 

Cti ~ Pi(l/(T - t,)) cos 2 (Zr s (r„)/ (r - r*)), 
Cei ~ m/(T - r*)) sin 2 {lr s {T,)l (r - r*)), 
Ctei ~ Wfo - n))sm(2lr s (n)/( To - r*)), (8) 

for adiabatic initial conditions. 

The oscillatory contributions to the C;s come only from the physics of the photon-baryon fluid, which is the same 
regardless of initial power spectrum. Therefore we may try to produce a non-oscillatory Cti spectrum with an 
oscillating primordial spectrum. We will show this in the next section. 

We see that Cti, Cei and Ctei all oscillate with the same frequency, but Cti, Cei are half a period out of phase 
and Ctei is a quarter of a period out of phase with both Cti and Cei- Therefore we cannot erase oscillations in more 
than one of these spectra at a time simply by changing the initial power spectrum of perturbations. If we produce 
a non-oscillatory spectrum of temperature perturbations, we must necessarily produce enhanced oscillations in the 
polarization and temperature polarization cross-correlation spectra. The question arises, however, of whether these 



enhanced oscillations would be observable. We look at this in section VI. 

The effect on the matter power spectrum is much simpler. As we look at this in Fourier space and we are only 
considering linear behaviour, the Fourier modes do not couple and we can relate the spectrum today directly to the 
primordial spectrum with a transfer function: 
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P f (k)=T a (k)P i {k), 



(9) 



where Pf(k) = \5f{k)\ 2 is the power spectrum today. 

Transfer functions for all species fall off for k > k eq as perturbations inside the horizon decay during the radiation 
dominated era. The baryon transfer function has oscillations inside the horizon before and during recombination, due 
to the oscillations in the photon-baryon fluid as described above, while the CDM transfer function is smooth since 
CDM interacts only gravitationally. After recombination however, the baryons stop interacting with the photons, 
effectively have no pressure, and only interact gravitationally, falling into the gravitational wells created by the CDM. 
The CDM and baryons therefore end up with the same transfer function today. In models with low baryon density 
the transfer function is smooth today, for models with high baryon density the oscillations may not have disappeared 
by the present day. In general, there could also be a hot or warm component to the dark matter, and these would also 
tend to damp the oscillations. For models in which the present day transfer function has no oscillations, oscillations 
in the primordial spectrum would still be present today. If the transfer function still has oscillations, eg due to a large 
baryon fraction, we would expect them to be at least partially cancelled by the primordial oscillations. In Section ^ 
we calculate the present day spectrum and discuss whether these oscillations are observable. 

We note that while we focus on the standard adiabatic picture here, one could construct a very similar discussion 
for the isocurvature case. Although specific aspects of the oscillations are different for isocurvature perturbations, the 
tendency for the fluid to oscillate and the resulting impact on the different observations is similar. In fact, a more 
general analysis would include a combination of all possible perturbation modes, as was done recently in ]l9| . It is 
clear from the results in [jl9f that the differentiating power of the polarization is present even when general mixtures 
of perturbation modes are considered. Introducing a wider fitting parameter space will naturally require more precise 
data to pin down the parameters, but we expect the main thrust of our analysis to persist under a more general 
multi-mode treatment. 



III. INITIAL POWER SPECTRUM 



There has been much interest recently in models of inflation which produce a broken-scale-invariant spectrum of 
perturbations. There is a model which produces oscillations in the primordial spectrum from an oscillating inflaton 
field [jl7 18 1, and there has been work on models with two fields in a potential with a step, producing oscillations in 
the spectrum below some scale |2(],|2l]]. Also multiple periods of inflation could produce oscillations in the spectrum 

We have chosen not to use a specific model to produce the initial spectrum, but to use a spectrum which completely 
cancels oscillations in the temperature spectrum created during the tight coupling era, as a worst case scenario, as 
this would be most confused with active models. The two models in this paper are examples for illustration of the 
ideas and analysis of how well they fit the data. They are intended to demonstrate some general features of this type 
of model. While Barrow and Liddle |23| argue that it might not be possible to construct a realistic inflaton that can 
produce the sort of primordial spectrum we discuss here, they only consider single-dimensional inflation space. If the 
data really do not show oscillations it is quite likely that some inflation enthusiast will find a model that tries to evade 
the secondary oscillation test. In any case our discussion is valid for the more general question passive vs. active 
perturbations. 

We have looked at two models: one with a cosmological constant and a low value of Ho (Model 1) and one with 
a higher Hubble constant and a tilted initial spectrum (Model 2). Model 1 has fib = 0.061, Q m = 0.32, fl\ = 0.68 
and h = 0.47. Model 2 has = 0.05, fl m = 0.3, ft\ — 0.7 and h — 0.7. Figure |l| shows the initial spectra with 
oscillations for the two models. Also shown are non-oscillatory spectra for comparison. 
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FIG. 1. Initial spectra with oscillations for Model 1 (left) and Model 2 (right). The dashed lines are the scale-invariant 
spectra. 



IV. TEMPERATURE ANISOTROPY SPECTRUM 



First we will compare the models to existing data on the microwave background. We have chosen to look only at 
high / since the spectra agree at low I. Figure |] shows the usual and modified Cti spectra for the two models, with all 
the data points currently available with an effective / of over 180. This includes Saskatoon fl24| , MSAM |2j|, Python 
HH], MAT H, CAT H, OVRO @, SUZIE @, Boomerang § and MAXIMA JToJ data. 

In order to compare the predictions from the models to experimental data we use the package Radpack 32 1 . This 
calculates the \ 2 between data and model in band powers, using a transformation of the C;s into approximately 
Gaussian distributed quantities [p3|. 




200 400 600 800 200 400 600 800 

1 1 

FIG. 2. Temperature anisotropy spectra for models with and without oscillations shown with current data with effective I 
of 180 or higher. Left: Model 1, right: Model 2 



Model 1 with a scale invariant initial spectrum has a % 2 per degree of freedom of 2.9 and the same model with an 
oscillating spectrum has \ 2 P er degree of freedom of 3.8. Model 2 has % 2 per degree of freedom of 0.94 with a non- 
oscillatory initial spectrum and 0.95 with an oscillating primordial spectrum, so fits the data very well. Introducing 
oscillations into the primordial power spectrum does not alter the level of fit appreciably. Clearly we cannot say 
anything yet about the presence or absence of oscillations. 

Future observations should be able to tell if there are oscillations in the CMB temperature spectrum. Figure |^ shows 
crrorbars (including the beam and uniform, uncorrelated noise, assuming foregrounds have been removed) predicted 
for the PLANCK satellite for Cti spectra with non-oscillatory spectra. Also shown are the spectra produced by an 
oscillatory primordial spectrum. 
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The errorbars are calculated in the standard way, assuming that the a\ m are Gaussian distributed and that there 
is spatially-uniform, Gaussian noise [j35|]36| l. The ler error bar in Cti is ACrz = \j2j (21 + l)f s k y (CTi + lo^B^ 2 ) 
where lotB 2 = ^2 c WT(c)Bfr c ): w t( c ) * s ^ ne var i ance of the noise in the temperature measurement in channel c, and 
B/( c ) ~ cxp (— l 2 cr 2 c \/2) is the Legendre transform of the beam for channel c. The angular resolution of the beam is 

given by ov c ) = ^/ t0 / lm ( c )/v / 81n2 where Of w hm is the full width at half max of the beam in radians. Incomplete sky 
coverage is partially taken into account by the factor f s k y which is the fraction of sky used to estimate the Cti ■ For 
PLANCK we have used = (O.OllfiK) 2 , combining the 143 GHz and 217 GHz channels at full width half max 
of 8.0' and 5.5' respectively, and f sky = 0.33 |§7). For the MAP satellite @ we used ujj, 1 = (0.106^if) 2 , with the 
highest frequency channels (40 GHz, 60 GHz and 90 GHz) at full width half max of 0.47°, 0.35° and 0.21° respectively 
and fsky = 0.66. The estimates of the errors are averaged in arbitrary bins to make the plots clearer. 

In order to determine whether or not future experiments will be able to distinguish between two models we need 
a different measure from the \ 2 statistic used to compare theory with observation. Since now we have two models 
which we wish to compare on an equal footing, we need a statistic that is symmetric in the two theories. We 
consider simulating Cs from each theory with mean C and standard deviation AC;. We use a Normal distribution 
to approximate the distribution of the Cs since we are only interested in the higher Is. 

— 1 -~-2 

There is a standard way to test if two samples have the same underlying mean: if C and Cj are both Gaussian- 
distributed with the same mean but possibly different variances then their difference will be Gaussian-distributed 
with mean zero and variance (AC/) 2 + (AC 2 ) 2 . So for each I we should use the statistic 

Y 2m = ( C ' ~ Cl I (10) 

X { ) (AC/) 2 + (AC 2 ) 2 [ ' 

and sum them over I. 

In fact we can use this idea to give a measure of how likely it is to confuse two models without performing simulations. 
The expectation of x 2 (0 is 

dl = 1+ (AC})*+ C (ACfr 

We will interpret the sum over the di in terms of a % 2 distribution, in order to find a probability with which the two 
theories could be confused. We define d to be the average of the di- Note that if the two theories are exactly the 
same, d = 1 as expected. We can use this measure to see how similar two theories are ideally, that is in the limit of 
no instrumental effects. To look at how well a particular experiment will distinguish the theories we put the noise, 
beam and fractional sky coverage effects into the spread in the C;s, the AC; as described above. 

Comparing the Model 1 C spectra with and without oscillations, using the MAP noise and beam parameters, gives 
a value for d of 2.6, which if interpreted as a x 2 per degree of freedom would mean that confusing the two models 
is ruled out at a a level of 43. For Model 2 d is 1.2, which corresponds to a a level of 3.9. The Model 2 spectra 
are harder to distinguish, since having a tilted initial spectrum suppresses the oscillations in the Cs. The PLANCK 
satellite will be more discriminating: with PLANCK characteristics d = 4.8 for Model 1 and 5.0 for Model 2. 
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Attempts have been made ]39|-|4l| to estimate the effect of foregrounds on our ability to measure the angular power 
spectrum. We can use the degradation factors for the temperature error bars in |fl[] to calculate d for the satellites 
taking into account their foreground models. The degradation factors Di increase the errorbars, ACf G = DiACi. We 
hope that the estimates for D\ given in |fjd|| would not vary too much with different underlying cosmological models. 
We use an upper limit D for the D; to give some room for variation. 

This gives us 

d FG = l + (d-l)/D 2 (12) 

For the main foreground scenario in D\ < 1.01 for all I values and both satellites which means that d is essentially 
unchanged. For their more pessimistic model Di < 2 which leads to d FG > 2.0 for PLANCK for both models, and 
d FG > 1.4 and d FG > 1.0 for MAP for Model 1 and Model 2 respectively. MAP will not be able to tell whether or 
not a tilted model has oscillations, but PLANCK can distinguish in both cases we consider here. 

To get the above values for di, we have not used I values below 50, as Ci are not Gaussian distributed for low I. 
The models we are interested in do not differ at those I, so d would decrease if we included them, but not enough to 
make a difference to the conclusions we draw. We have calculated each di separately, without binning, but this would 
make errorbars smaller and therefore d even bigger. 



V. MATTER POWER SPECTRUM 



In this section we look at how the altered primordial spectrum affects the present day matter power spectrum and 
whether these changes are observable. We would like to compare the model predictions for the matter power spectrum 
with current data, to see what we can already say about models with oscillations in the primordial spectrum. In the 
currently favoured models the transfer functions are smooth, and any features in the initial spectrum would show up 
in the present day matter power spectrum. For such models the power spectrum of the matter density field is the 
best probe of possible oscillations in the initial spectrum. 




k (h Mpc-') k (h Mpc- 1 ) 

FIG. 4. Present day matter power spectrum from the models with oscillations in the primordial spectrum, shown with data 
from Q. Left: Model 1, right: Model 2. 

Figure |4| shows the matter power spectrum from the models with and without oscillations and data from p2] ] . We 
use the decorrelated linear spectrum which has effectively uncorrelated points. It is important to have independent 
points since oscillations would be smoothed out in correlated data. We have calculated the x 2 f° r our models binned 
in the same way as the data. For Model 1 with a scale invariant initial spectrum the y 2 per degree of freedom is 3.2 
and with oscillations in the initial power spectrum it is 3.4. For Model 2 with oscillations the \ 2 P er degree of freedom 
is 1.2 with and without oscillations. This model fits the data well but we cannot say whether or not the spectrum 
has oscillations. 

We can look at averaged values, such as cr 8 . For Model 1 this is 0.53 and 0.56 without and with oscillations in the 
primordial spectrum, normalised to COBE. For the Model 2 the values are 0.68 and 0.74 without and with oscillations. 
We see that oscillations in the spectrum have very little effect on the value of as- We would expect this, as we have 
not changed the overall tilt of the spectrum. 

It does not seem that the current data on the matter power spectrum can tell us much about the presence or 
absence of oscillations as we have discussed here. In order to draw conclusions about this we need data points with 
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a smaller coherence length. We hope that new data from redshift surveys will constrain the power spectrum much 
more than has been possible so far E| , 

One possible complication would be if the preferred model turns out to have an oscillatory matter transfer function 
(eg because the matter has a high proportion of baryons). Even in that case, with good enough data it might be 
possible to tease apart the different effects. 



VI. POLARIZATION IN THE MICROWAVE BACKGROUND 

As we saw in Section [S| the power in the polarization anisotropy oscillates out of phase with that of the temperature 
anisotropy, so the extra power on certain scales in the initial spectrum leads to enhanced polarization oscillations. 
Here we show examples for our models. Note that active models have completely different polarization spectra so 
there is no opportunity for confusion [^j4^| . 

Figure |5| shows the angular power spectra for polarization for the two models using a scale invariant primordial 
spectrum with la errors Jll| expected from the PLANCK satellite. The parameters used are the same as for the 
temperature spectra earlier. 




FIG. 5. Angular polarization power spectra for models with and without oscillations, with errorbars forecast for the PLANCK 
satellite. The errorbars are shown for the models with scale-invariant initial spectra. Left: Model 1, right: Model 2 

The figures show clearly the enhancement of the polarization spectrum. To see if these differences are observable 
we have used the same measure d as in Section ^ to compare two models. Using errorbars predicted for PLANCK 
d — 6.5 for Model 1 and 2.8 for Model 2. For MAP the values are 1.01 and 1.00 respectively. So we would not consider 
these two curves different under MAP observations, but PLANCK should be able to distinguish them. 

In the main (supposedly most realistic) model of polarization foregrounds from jil| the degradation factor D\ for 
errorbars on the angular polarization spectrum is less than 2 for PLANCK at all I, which gives d FG > 2.4 for Model 
1 and d FG > 1.5 for Model 2. So PLANCK will be able to tell these models apart even in the presence of this level of 
foregrounds. The more pessimistic model gives Di < 10 for PLANCK at I > 50, giving d FG > 1.06 for Model 1 and 
d FG > 1.02 for Model 2, so if the foregrounds are this bad even PLANCK will not separate these spectra. 

Figure ^ shows the temperature-polarization cross-correlation for the various models. Here again we can see the 
shift in oscillations produced by altering the primordial spectrum. Looking at the cross-correlation power should be 
more valuable than the polarization on its own, for two reasons. The first is that the noise begins to dominate at 
higher I for temperature than for polarization, so we will get a clearer signal for the cross-correlation than for the 
polarization alone. Using errorbars predicted for PLANCK, d = 5.0 for Model 1 and 3.2 for Model 2. For MAP 
the values are 1.09 and 1.01 respectively. These are too low to allow MAP to distinguish between the models, but 
PLANCK would not be confused. We see that for these models the cross-correlation is in fact not such a good 
discriminator between models with and without oscillations in the temperature spectrum. This is because although 
the errorbars are smaller, the models are closer together. The effect of oscillations in the primordial spectrum is clearer 
in the present day polarization spectrum than in the temperature-polarization correlation since in the polarization 
the oscillations are completely out of phase, but in the cross-correlation they are only half a period out of phase. 

There are no estimates of the degradation factor for the cross-correlation in |ilj . We make a rough estimate using 
the fractional errors on C/s with foregrounds obtained by p9| . These give degradation factors slightly larger than 
the main foreground scenarios in j4l[] for the temperature and polarization spectra as found be PLANCK. We find 
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d FG > 1.6 for Model 1 and d FG > 1.4 for Model 2, hence we estimate that with this level of foreground contamination 
the PLANCK data should discriminate between these models with and without oscillations. 




i i 
FIG. 6. Temperature-polarization cross-correlation for models with and without oscillations, with errorbars forecast for the 
PLANCK satellite. The errorbars are shown for the models with scale-invariant initial spectra. Left: Model 1, right: Model 2 

The second reason for looking at the cross-correlation power is that a strong prediction for active models is that due 
to their incoherence the cross-correlation in models where active sources play the major role in structure formation 
will be close to zero. We have calculated d for the cross-correlation using the models with oscillations in the primordial 
spectrum and a 'model' where all the Ctei are zero, to see if the cross-correlation could be consistent with zero. The 
value of d for Model 1 is 15.4, and for Model 2 the value is 13.5. We would not consider either of these cross-correlation 
spectra to be consistent with a vanishing spectrum. Therefore an observation of a smooth temperature anisotropy 
spectrum and vanishing cross-correlation would be strong evidence against passive fluctuations. 



VII. DISCUSSION 



We have investigated the secondary oscillation test, which is recognized as the most fundamental way in which 
scenarios with perturbations from inflation can be falsified. While current data are not good enough for the test 
to be applied, we look toward the future when much better data will be available. We have studied particular 
situations where the inflationary model conspires to evade the test by building oscillations into the inflaton potential 
to effectively 'cancel out' the natural oscillations that would appear in the CMB anisotropics. We have shown that 
such conspiratorial models will produce the characteristic oscillations in other observable spectra, and thus will be 
unable to evade the test. 

One might think that reionization could smooth the peaks of the angular spectra and confuse us, but reionization 
also produces a different signal in the polarization. Reionization damps the temperature perturbations on all scales 
smaller than the sound horizon at the start of the reionization. Homogeneous reionization causes approximately 
exponential damping in I and so cannot alter the oscillations in the spectrum. In order to do this one would have 
to have damping starting on smaller scales than the sound horizon at recombination; clearly the reionization scale 
is larger than that of recombination. Reionization damps the polarization spectrum as well but also produces more 
polarization on large scales due to rescattering at the second 'last scattering surface', resulting in a boost at low I in 
the polarization angular spectrum. The oscillations in the spectra might be washed out enough to be unobservable 
if the optical depth to reionization is large enough, but increasing optical depth increases the extra power in the 
polarization at low I [^| , providing a strong indication of the presence of reionization. So inflationary perturbations 
with homogeneous reionization could not evade the secondary oscillation test. 

The situation might be different in non-standard reionization where the visibility function is broader, but this 
causes a large amount of damping and is probably ruled out already by the temperature observations. In inhomo- 
geneous reionization models the first order temperature spectrum comes entirely from the homogeneous reionization 
background. The second order effect due to the inhomogeneities is too small compared to the first order to make any 
difference to the peaks (Hill ■ 

Work has been done recently on models with both active and passive fluctuations p9| ]. The application of the fine 
tuning we have done here to those models represents an even more exotic corner of parameter space which we have 
not examined. 







The matter power spectrum is affected by a number of issues (such as the baryon fraction) which make the inter- 
pretation of the presence or absence of oscillations less than clear-cut. However, we have shown that the polarization 
power spectrum and the temperature-polarization cross-correlation will offer a much more straightforward test, as- 
suming the cosmic signal can be separated from the foregrounds. It is these spectra we are relying on to draw our 
optimistic conclusions. 

There is no question that the models we do consider are viewed by most cosmologists (including ourselves) as 
extremely unnatural and fine-tuned. A very reasonable stand might be that good taste alone would require us to 
abandon an inflationary origin to the perturbations rather than embrace the conspiratorial models we study here. 
The main point of this work is that modern datasets will leave nothing to good taste. By the time the new data is 
in, even the exotic part of model space we explore here will be exposed to confrontation with observations. 
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